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Abstract -Recent experimental and theoretical studies focus on spin-mediated heat currents 
at interfaces between normal metals and magnetic insulators. We resolve conflicting estimates 
for the order of magnitude of the spin transfer torque by first-principles calculations. The spin 
mixing conductance G 1 ^ of the interface between silver and the insulating ferrimagnet Yttrium 
Iron Garnet (YIG) is dominated by its real part and of the order of 10 14 f2 _1 m -2 , i.e. close to 
the value for intermetallic interface, which can be explained by a local spin model. 
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Introduction. — It has recently been reported that 
the magnetism of insulators can be actuated electrically 
and thermally by normal metal contacts Q3H]- The mate- 
rial of choice is the ferrimagnet YsFesO^ (YIG), because 
of its extremely small magnetic damping [3H3]. The low- 
lying excitations of magnetic insulators are spin waves, 
which carry heat and angular momentum [B]. Existing 
experiments use Pt contacts, which by means of the in- 
verse spin Hall effect are effective spin current detectors 
[7]. Slonczewski [5] reports that the thermal spin transfer 
torque in magnetic nanopillars [5|10| can be much more ef- 
ficient than the electrically generated spin torque in metal- 
lic structures. 

The electrical and thermal injection of spin and heat 
currents into insulating magnets is governed by the spin 
transfer torque at the metal | insulator interface [TT1IT2] . 
which is parameterized by the spin-mixing conductance 
G^ = e 2 Tr (i - rfrj /h, where I and r CT are the unit 
matrix and the matrix of interface reflection coefficients 
for spin a spanned by the scattering channels at the Fermi 
energy of the metal [T3]. Crude approximations such as 
a Stoner model with spin-split conduction bands |TT] and 
parameterized exchange between the itinerant metal elec- 
trons and local moments of the ferromagnet [T1[8iri2] have 
been used to estimate G^ for YIG interfaces^. Experi- 



ments and initial theoretical estimates found very small 
spin torques that arc at odds with Slonczcwski's predic- 
tions [8]. 

Here we report calculations of the spin mixing conduc- 
tance for the Ag|YIG interface based on realistic electronic 
structures. Silver is a promising material |14j for non-local 
spin current detection |15j , which should be more efficient 
than the inverse spin Hall effect in nanostructures. We 
demonstrate that the calculated G^ for the Ag|YIG inter- 
face is much larger than expected from the Stoner model 
and better described by local-moment exchange fields. 

Free-electron model. — We start with a reference 
structure consisting of an Ag|FI| Ag(OOl) junction in which 
the ferromagnetic insulator (FI) is modeled by a spin- 
split vacuum barrier, i.e., the free-electron Stoner model. 
The vacuum potential is chosen to be spin-split by 0.3 
and 3.0 eV, whereas the barrier height is adjusted to 0.3, 
1.4, 2.6 and 2.85 eV, respectively. The barrier thick- 
ness (1.2 nm) is chosen here such that electron transmis- 
sion is negligible. Table [1] lists the corresponding G^ of 
Ag|FI| Ag. Both Re G^ and Im G^ decrease with increas- 
ing barrier height, as expected |11| . 

Band structure. — We calculate the electronic struc- 
ture of YIG using the tight-binding lincar-muffin-tin- 



' a ' E-mail: kexia@bnu.edu.cn 
lr The spin mixing conductance is governed by the reflection coeffi- 
cients only and remains finite when the transmission coefficients van- 



ish. This is not a breach of the scattering theory of transport, since 
the incoming and outgoing scattering states are well defined as prop- 
agating states in the metallic contacts. 
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Table 1: Spin-dependent and spin mixing conductances of a 
Ag|FI| Ag(OOl) junction with different barrier heights and spin 
splitting A = 0.3 and 3.0 eV. The mixing conductances for 
the (111) orientation differs by less than 20%. The Sharvin 
conductance (G Sh ) of Ag(001) is 4.5 x 10 14 m~ 2 . 



Barrier G r /G sh 


G^/G sh 


RcG n /G sh 


ImG^/G' 


A = 0.3 eV 








0.3 6.3E-5 


5.1E-6 


0.009 


-1.1E-1 


1.4 3.3E-8 


7.1E-9 


0.003 


-7.4E-2 


2.6 3.5E-9 


1.3E-9 


0.001 


-4.0E-2 


2.85 0* 





0.001 


-5.1E-2 


A = 3.0 eV 








0.3 7.0E-6 





0.15 


-0.45 


1.4 7.4E-10 





0.08 


-0.35 


2.6 





0.05 


-0.28 


2.85 





0.04 


-0.27 


* means a transmission probability of less than 10 



orbital code in the augmented spherical wave approxi- 
mation as implemented in the Stuttgart code j!6H18j us- 
ing the generalized gradient correction (GGA) to the lo- 
cal density approximation (LDA). The cubic lattice con- 
stant a = 12.2 A is chosen 1.6% smaller than the exper- 
imental one [T5]. We use 136 additional empty spheres 
(ES) for better space filling and reduced overlap between 
neighboring atomic spheres. YIG is a fcrrimagnetic insu- 
lator with band gap of 2.85 eV [5D1[3T]. Magnetism is 
carried by majority and minority spin Fe atoms (tetrago- 
nal Fe(T) and octahedral Fe(O) sites in FigOJa), respec- 
tive.) with a net magnetic moment of 5/xb per formula 
unit [T9"l |2"2"H2~4"] . The magnetic moments are 3.95 and 
—4.06 fiB for majority and minority spin Fe atoms, respec- 
tively. Both Y and O atoms show small positive magnetic 
moments of 0.03 and 0.09 fiB, respectively, while those on 
the empty spheres do not exceed 0.007 /is. The common 
problem of density-functional theory to predict the energy 
gap of insulators can be handled by an on-site Coulomb 
correction (LDA/GGA+U) [25J.26 or a scissor operator 
(LDA/GGA+C) [27]. Figure []]» is a plot of the band 
structure of GGA with a fundamental band gap of 0.33 eV 
between the valence band edge of the majority-spin chan- 
nel and conductance band edge of minority-spin channel. 
The GGA+C method can be used to increase the band 
gap depending on the scissor parameters C. A GGA+C 
band structure with a band gap of ^1.25 eV is shown in 
Fig. [He). The GGA+U method applied to the YIG band 
structure using the parameters from Ref. [25[|2l)] leads to 
the band structure plotted in Hid) with the same energy 
gap - 1.25 eV. 

While a visual comparison of the band structures in 
Figs. [Uc) and (d) assures the equivalence of the two meth- 
ods, we can assess the differences quantitatively by com- 
paring the effective masses at the band edges as shown in 



Table 2: Band gap (E g ) and effective masses (in unit of m e ) of 
the band structure of YIG at the V point as calculated by the 
GGA, GGA+U, and GGA+C methods. CB and VB denote 
conductance and valence bands, respectively. 





E g (eV) 


Majority-spin 


Minority-spin 






VB 


CB 


VB 


CB 


GGA 


0.33 


0.10 


0.52 


0.40 


0.17 


GGA+U a 


1.25 


0.13 


0.60 


0.37 


0.19 


GGA+C b 


1.25 


0.17 


1.00 


0.31 


0.27 


GGA+C C 


1.4 


0.17 


1.00 


0.28 


0.31 


GGA+C d 


1.4 


0.18 


1.46 


0.25 


0.25 


a U = 3.5 eV, J = 0.? 


; eV 








b C{Fe,Y) 


= 6.1 eV,C(ES) 


= 3.05 eV 






c O(Fe,Y) 


= 7.2 eV 










d C(Fe,Y) 


= 7.5 eV,C(ES) 


= 3.75 eV 







Tablc[2l For band gaps of ^ 1.25 eV the effective mass at 
the conductance band edge of majority-spin as obtained 
by the GGA+U and GGA+C methods differ by up to 
67%. This seems significant, but the effects on the mixing 
conductance, which is the quantity of our main interest 
here, is small, as discussed in the next section. 

Ag|YIG interface. — Wc study the spin mixing con- 
ductance in Ag|YIG|Ag with a 3x3 and 6x6 lateral su- 
percell of fee Ag to match a cubic YIG unit cell along the 
(001) and (111) directions with lattice mismatch ~ 1%. 
Interfaces can be classified according to their magnetic sur- 
face properties into three types of terminations. For the 
(001) texture, one cut is terminated by Y as well as ma- 
jority and minority spin Fe atoms with compensated mag- 
netic moment ( ' YFc-termination" ) . Another cut yields 
only majority Fe atoms at the interface ( "Fe-tcrmination" ) 
with total magnetic moment of 7.90 /j,b per lateral unit 
cell. The third interface covered by O atoms is obtained 
by removing Fe and Y atoms from the YFe-termination. 
The oxygen layer is separated from adjacent Fe atoms by 
only ~ 0.3 A. Including the latter, the "O-termination" 
also corresponds to a net interface magnetic moment of 
7.90 /J.B- The interfaces for the (111) direction can be clas- 
sified analogously. The "YFe-termination" cut has now 
a net interface magnetic moment of 23.70 /is. The Fe- 
tcrmination contains now minority-spin Fe atoms with net 
magnetic moment of — 16.24/is, while the O-terminated 
surface has the same magnetic moment when including 
the shallowly buried Fe layer. 

We chose a YIG film of 4 unit cell layers, because its 
electric conductance does not exceed 10~ 10 e 2 /h per unit 
cell. G^ is therefore governed solely by the single Ag|YIG 
interface. 

First, we inspect G^ of Ag|YIG interfaces computed 
with and without scissor corrections. We find that the 
difference of ReG^ is as small as 21% when increasing 
the band gap of YIG from its GGA value of 0.33 eV to 
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( c ) GGA+C ( d ) GGA+U 



Majority-spin ^ Minority-spin Majority-spin Miiority-spin 




Fig. 1: (a): 1/8 of the cubic YIG cell; the full structure can be 
obtained by symmetry operations. Here, Fe(T) and Fe(O) are 
Fe atoms at tetragonal and octahedronal sites, respectively, (b- 
d): Band structures of YIG with GGA, GGA+C, and GGA+U 
method with band gap of 0.33, 1.25, and 1.25 eV, respectively. 

a GGA+C band gap of 2.1 cV as shown in Table [3] We 
conclude that the precise band gap is a parameter that 
hardly affects the G^ of the Ag|YIG interface. 

Table 3: Spin mixing conductance of Ag|YIG(001) with differ- 
ent YIG band gaps modulated by the GGA+C methods. We 
pin the Fermi level of Ag at mod-gap of YIG. 0.33 eV is the 
band gap of GGA (without a scissor operator). 



E g (cV) 


0.33 


0.65 


0.95 


1.4 


1.8 2.1 


RcG^ 














3.46 


3.94 


3.43 


3.01 


2.82 2.74 



By scanning the Fermi energy of Ag (or the YIG work 
function), we can obtain information similar to that when 
varying the band gap. Scanning the Ag Fermi energy from 
the valence to conductance band edges for a band gap of 
1.4 eV, we obtain results equivalent to a mid-gap Fermi 
energy and band gaps varying from zero to 2.8 eV, but 
without changing the details of the band dispersion. In 
Figure [5]we plot the mixing conductance of Ag|YIG(001) 
with YFe termination as a function of YIG's work func- 
tion. Here we consider two kinds of band dispersions with 
the same band gap of 1.4 eV obtained by different scissor 
operator implementations as shown in TablcJSJ We find 
that the mixing conductance does not depend sensitively 
on (i) the YIG work function or interface potential bar- 
rier as well as (ii) the band dispersion when fixing the 
Fermi energy of Ag in the middle of the band gap of YIG; 
the difference in effective mass of 46 % causes changes in 




_1 i I i I i I i I i I i i L_ 

0.2 0.4 0.6 8 1.0 1.2 1.4 0.0 

Energy with respect to YIG (eV) CB VB 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 0.0 

Energy with respect to YIG (eV) CB VB 



Fig. 2: Effect of band dispersion and band alignment on 
the spin mixing conductance of Ag|YIG| Ag(001) with YFe- 
termination. We use YIG with same band gap of 1.4eV with 
different implementations of scissor operator (a) C(Fe,Y) = 
7.2 eV; (b) C(Fe,Y) = 7.5 eV, and C(ES) = 3.75 eV to see the 
effect of band dispersion. We fix the Fermi energy of Ag while 
scanning the YIG work function. 

RcG^ of only 13%. These deviations are within the error 
bars due to other approximations (see below). We there- 
fore conclude that the transport properties in the present 
system are sufficiently well represented by the scissor op- 
erator or on-site Coulomb correction methods for the gap 
problem. 

Besides the band alignment discussed in the previous 
paragraph, two more properties are difficult to compute 
self-consistcntly for large unit cells, viz. the atomic inter- 
face configuration and the ferromagnetic proximity effect: 

(i) : We determine the distance between Ag|YIG by mini- 
mizing ASA overlap while keeping the space filled. We es- 
timate that the differences in G^ for configurations with 
maximum and minimum ASA overlap is less than 30% 

(ii) : We assess the ferromagnetic proximity effect by us- 
ing the self-consistent electronic structure of Ag atom at 
the Ag|Fe interface. We find that the Ag atoms closest 
to Fe acquire a magnetic moment of 0.025 \xb and the ef- 
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Fig. 3: Spin mixing conductance of Ag|YIG(001) (left) and 
Ag|YIG(lll) (right) with a YIG band gap of 1.4 eV. We fix 
the Fermi energy of Ag while scanning the YIG work function. 

feet is observable up to the 4th Ag layer. The spin mixing 
conductance is found to be enhanced by about 10% in this 
system. In the following we disregard such an effect. From 
various checks of these and other issues, the magnitude of 
a possible systematic error in the mixing conductance is 
estimated to be < 40%. 

Results. — Fig. [3] summarizes our results for G^ 
of Ag|YIG(lll) and Ag|YIG(001) junction with different 
YIG interface-terminations. We find G^ ~ 10 14 f2 _1 m~ 2 
for both Ag|YIG(lll) and Ag|YIG(001), with the real 
part dominating over the imaginary one, in stark con- 
trast to the Stoncr model (cf. Table 1). G^ depends 
only weakly on exposing different YIG(lll) surface cuts 
to Ag, which we attribute to the homogeneous distribution 
of magnetic atoms. The YFe termination of YIG(OOl) is a 
nearly compensated magnetic interface, but we still calcu- 
late a large spin mixing conductance. Finally, our results 
are two orders of magnitude larger than the experimental 
value found for Pt|YIG(lll) [T]! The difference between 
Ag and Pt cannot account for this discrepancy: one would 
rather expect a larger G^ for Pt because of its higher con- 
duction electron density. 

The difference between the Stoner model and the 
first-principles calculations indicate that the spin-transfer 
torque physics at normal metal interfaces with YIG is very 
different from those with transition metals. Spin-transfer 
is equivalent to the absorption of a spin current at an in- 
terface that is polarized transversely to the magnetization 
direction. Magnetism in insulators is usually described in 
a local moment model. The physical picture of spin trans- 
fer appropriate for metals, viz. the destructive interference 
of processing spins in the ferromagnet, then obviously fails. 
When the spin transfer acts locally on the magnetic ions, 
we expect no difference for the spin absorbed by a fully 
ordered interface with a large net magnetic moment or a 



Fig. 4: ReG n as a function of interface magnetic moment den- 
sity in Ag|YIG compared with that of Fe atoms at the Ag|Vac 
interface. Insert: Crystal-plane resolved spin density {cr y ) in ar- 
bitrary units for Ag|YIG(001) with YFe-termination when fully 
polarized electrons are injected from the Ag side at mid-gap en- 
ergy. The maximum magnetic moment density is 39/is/nm 2 as 
estimated from a full monolayer of Fe at the interface, in which 
the Fe atoms adopt the Ag structure with magnetic moment 
of 2.81/iB- 

compensated one, in which the local moments point in op- 
posite directions, as is indeed born out of our calculations. 

In order to test the local moment paradigm, we consider 
non-conducting Ag|Vac(4L)|Ag(lll) junctions. We now 
sprinkle one vacuum interface randomly with Fe atoms. 
At low densities the Fe atoms are weakly coupled and 
form local moments. The electronic structure is gener- 
ated using the Coherent Potential Approximation (CPA) 
for interface disorder. A 10 x 10 lateral supercell with 
100 atoms in one principle layer is used to model a mag- 
netic impurity range from 1% to 80%. The high den- 
sity limit is a monolayer of Fe atoms in the fee struc- 
ture: Ag|Fe(lL)|Vac(3L)|Ag(lll) with total magnetic mo- 
ment of 2.81 fiB per Fe atom. So, the maximum mag- 
netic moment density here is 39 fiB / nm 2 . The results for 
the mixing conductances is summarized in Fig. 01 We 
find that the ratio of G^ to the (Ag) Sharvin conduc- 
tances monotonically increases with the Fe density at the 
Ag|Vac interface. The increase is linear at small densi- 
ties and saturates around 30 /j,b/ nm 2 due to interactions 
between neighboring moments. We find that ReG^ of 
Ag|YIG and Ag|Fe|vacuum agrees well for corresponding 
Fe densities at the interface, in strong support of the local 
moment model. 

Since the mixing conductance is dominated by the local 
moments at the interface, we understand that the results 
are relatively stable against the difficulties density func- 
tional theory has for insulators. The variation of the band 
gap of the insulator as well as the band alignment with 
respect to normal metal changes the penetration of the 
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spin accumulation, but since only the uppermost layers 
contribute this is of little consequence. 



Table 4: G n of a disordered Ag|YIG(001) interface with YFe- 
termination. Directional disorder is introduced by flipping 
three majority Fe spins in the 2x2 super cell. 





ReG^ (10 14 fi _1 m- 2 ) 


lmG n (lO^fi^m 


clean 


3.010 


0.302 


disorder 


3.145 


0.382 



Table |4] shows the effect of directional disorder 
of magnetic moments on the mixing conductance for 
Ag|YlG(001) with YFc-termination, for which the in- 
tegrated surface magnetic moment density is close to 
zero. Here, we use a 2 x 2 lateral Y1G supercell in 
which three magnetic moments arc flipped to a negative 
value, amounting to a total surface magnetic moment of 
— 23.7 Hb per lateral unit cell. The directional disorder 
of magnetic moments at the interface slightly enhances 
ReG^ (around 5%), as indeed expected from the local 
moment picture. 



the normal metal and YIG should enhance the spin mixing 
conductance. In Fig. [5] (a), we show that inserting Fe 
atomic layers indeed increases ReG^" by 40-65% up to 
the intermetallic Ag|Fe value, which is close to the Ag 
Sharvin conductance. 

In these calculation, the Ag reservoirs has been assumed 
to be ballistic. When the spin mixing conductance is small 
2 .relative to the Sharvin conductance, this is a valid approx- 
imation, but otherwise the diffusive nature of transport 
may not be be neglected. Since ReG^ turns out to be 
of the same order as G^ we have to introduce the diffu- 
sive transport correction as introduced by Schep et al. as 

i i i 



^E 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

VB Energy with respect to YIG (eV) CB 



Fig. 5: Spin mixing conductance of Ag|Fe„|YIG|Fe„|Ag(001) 
with YFe termination for (a) ballistic and (b) diffusive trans- 
port (i.e. in the presence of Schep correction), where n is the 
number for Fe monolayers inserted between Ag and YIG. 



Gti Gt-i- 



(1) 



The results are shown in Fig. E][b). We observe that the 
" Schep" correction enhances the spin mixing conductance 
by 20% for for the and about 90% for 4 nonolayers Fe 
insertions between Ag and YIG. 

The spin transfer can be maximized by a high den- 
sity of magnetic ions at the interfaces. In YIG we could 
not identify interface directions or cuts that are espe- 
cially promising, but this could be different for other 
magnetic insulator, such as ferrites [5]. Slonczewski [5] 
uses a local moment model with a somewhat smaller ex- 
change splitting (0.5 eV) than found here; when defined 

as A (r) = f nws (Vj c (i?,r) - Vl (i?,r)) p (R,r) dr, 

where p (r, r) is the density of the evanescent wave func- 
tion in YIG at mid-gap energy disregarding its spin split- 
ting [30], flws the Wigncr-Scitz sphere at the lattice site 
i?, and Vxc denotes the exchange-correlation potentials 
for spin-up (down) electronsthe exchange splitting felt by 
the Ag conduction electrons at the YIG interface is up to 
~ 3.0 eV. Since Slonzcewski focusses on the magnetiza- 
tion dynamics of the magnetic insulator we cannot carry 
out a quantitative comparison with his model here. 

Conclusion. — In conclusion, we computed the spin 
mixing conductance G^ of the interface between sil- 
ver and the insulating ferrimagnet Yttrium Iron Garnet 
(YIG). Re G n is found to be of the order of 10 14 fT 1 m" 2 , 
which is much larger than expected for a Stoncr model, 
which indicates the importance of the local magnetic ex- 
change field at the interface. On the other hand, G^ is 
not very sensitive to crystal orientation and interface cut. 
Re G^ can be enhanced to around 40-65% of the fully 
metallic limit by inserting a monolayers of iron between 
Ag and YIG. The discrepancy between the measured and 
calculated mixing conductance might indicate previously 
unidentified interface contaminations that, when removed, 
would greatly improve the usefulness of magnetic insula- 
tors in spintronics. 



Inserting a thin ferromagnetic metallic layer between We would like to thank Burkard Hillebrands, Eiji 
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Additional remark: After first submission of our 
manuscript arXiv:1103.3764, a manuscript was submitted 
and accepted by Physical Review Letters (Heinrich B. et 
al., Phys. Rev. Lett., 107 (2011) 066604) that reports a 
mixing conductance that is clearly enhanced compared to 
ref. [1], but still an order of magnitude smaller than our 
predictions. 
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